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The propagation of waves in highly inhomogeneous media is a problem of interest in multiple 
fields including seismology, acoustics and electromagnetism. It is also relevant for technological 
applications such as the design of sound absorbing materials or the fabrication of optically devices 
for multi-wavelength operation. A paradigmatic example of a highly inhomogeneous media is one in 
which the density or stiffness has fractal or multifractal properties. We investigate wave propagation 
in one dimensional media with these features. We have found that, for weak disorder, localization 
effects do not arrest wave propagation provided that the box fractal dimension D of the density 
profile is D < 3/2. This result holds for both fractal and multifractal media providing thus a simple 
universal characterization for the existence of localization in these systems. Moreover we show that 
our model verifies the scaling theory of localization and discuss practical applications of our results. 

PACS numbers: 05.45.Df.73.20.Fz, 05.60.Cd,46.40.Cd,46.65.+g 



Wave propagation [l|, y] in disordered or highly in- 
homogeneous media is a recurrent topic of research in 
physics. This is not surprising due to its broad range 
of potential apphcations. The propagation of electro- 
magnetic waves in highly inhomogeneous fractal media 
[S-Q was studied experimentally in [J, Q ■ Anderson lo- 
calization effects in acoustic waves have been observed 
in experiments |7| and numerical simulations [8|. Seis- 
mic waves [9| also propagate in a highly inhomogeneous 
medium as there is a growing consensus [10, \lM that the 
distribution of fractures and densities in the earth inner 
structure has a fractal distribution. A detailed under- 
standing of this problem is also relevant for the fabrica- 
tion of optically active devices for multiwavelength oper- 
ation [3|, the fabrication of a multiple gap and multiple 
pass band micro-strip resonator filter ^Q], the design of 
sound absorbing materials and the nondestructive char- 
acterization of fractured materials [12] • 

Rigorous theoretical results are known only for one- 
dimensional (Id) matter waves in correlated random po- 
tentials V{n). Kotani [13] demonstrated that Anderson 
localization - stop of diffusion and exponential decay 
of eigenfunctions - occurs for any energy and amount 
of disorder [l|, [^ provided that correlations decay as 
{V{n)V{0)) oc 1/n" with a > or faster. Anderson 
localization effects are also suppressed as the degree of 
differentiability of the potential increases 1J| . For other 
recent studies of localization in correlated Id potentials 
we refer to [la| . It is unclear to what extent these results 
still hold in the case of classical waves in a fractal-like 
media. 

The main aim of this paper is to answer this question. 
We investigate what are the most general circumstances 
in which a perturbation in a highly heterogeneous Id 
medium, described by a fractal or multifractal density, 
propagates without ever experiencing Anderson localiza- 
tion. 
Our main result is that a perturbation will propagate in- 



definitely in Id provided that the box fractal dimension 
of the medium density/stiffness fluctuations is D < 3/2. 
These results are still valid if the density profile is mul- 
tifractal. 
Monofractal density.- We investigate a discretized ver- 




FIG. 1. Average length {L{k)) of the medium over an interval 
of width k. By definition the slope of the curve is the box 
dimension D. The results correspond to a density profile given 
by a fBm with Hurst exponent H. H ranges from 0.1 to 0.8, 
in steps of 0.1, from botton to top. In all cases D = 2 — H. 
Inset: Generalized Hurst exponent h{q) (/i(2) = H) [231 that 
describes the scaling of a multifractal density profile generated 
by a multiplicative random cascade process (see text). 

sion of the Id scalar wave equation, 

p„(V'«+l - V'n) - Pn-l(V^n - V'n-l) + '^^V'ri = 0, (1) 

where a; are the eigenfrequencies of the system, /?„ = 
po(l + 'yn) the ratio between the stiffness and the density 
of the material and po is the average value. The fluc- 
tuations rjn are described by a fractal function with box 
fractal dimension D. For the sake of simplicity we study 



localization properties in a narrow region of frequencies 
outside the origin w w as a function of the box dimen- 
sion D and the system size L. We will restrict ourselves 
to the region of weak disorder. This a natural choice as 
our main motivation is to put forward a characterization 
for the existence of delocalized eigenmodes. As examples 
of density profiles with a monofractal spatial -labelled by 
n- distribution we employ: (i) the Weierstrass curve 
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= X]cos(7'=n/L + 0fe)/7(2-^)^ 



(2) 



with 7 > 1, 0fc is a box distributed random number 
e [0, 2tt] and n is a spatial variable. The box fractal 
dimension is D'^'^ = D ^Q,UM^ (ii) ^ fractional Brown- 
ian motion (fBm). This is a generalization of the stan- 
dard Brownian motion characterized by stationary incre- 
ments ~ though the process itself is non-stationary-, self- 
similarity, and a variance cr^ oc n^^ where H is the Hurst 
exponent. We note that, due to its self-similar character, 
every realization of the process is a fractal curve with 
box-dimension D = 2 — H [20|. In both cases in order 
to have a well defined continuous limit we rescale the po- 
tential such that the average is < 77 >= 0. The standard 
deviation ar^ = •\/< rf > controls the strength of disor- 
der. The box-dimension D for agiven H is obtained by 
the method presented in Ref. [l3|- It basically consists 
in calculating the average curve length {L{k)) over an 
interval of width k. If {L{k)) oc fc~^, then the curve is 
fractal with a box-dimension D. In Fig. [l]we show that 
for the fBm I?fBm = 2 — H. We compute the eigenmodes 
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FIG. 2. (color online) var(s), as a function of D for u ~ 0.1, 
(7,, =0.1 and different sizes L obtained by numerical diag- 
onalization of the wave equation Eq.([T]) for a density profile 
described by the Weiertrass function Eq.(l2]), and a fractional 
Brownian motion (see text for details). In both cases var(s) 
increases with the system size - a signature of Anderson local- 
ization - only for D > 3/2. This indicates that D = Dc = 3/2 
is the maximum box dimension for which a band of extended 
modes can occur. 

cj of Eq. Il]) by using standard numerical diagonalization 
techniques. For a given disorder a^i and frequency win- 
dow the number of eigenmodes obtained is at least 2x 10^. 



In order to determine the importance of Anderson local- 
ization effects we carry out a finite size scaling analysis of 
the spectrum [I^]. The variance var(s) of the eigenmode 
spacing distribution P{s) is chosen as the scaling variable 
[IS] . P{s) is the probability of finding two neighbouring 
eigenmodes at a distance s = (cj^+i — Wi)/A and 



var(s) = (s^) - (s)^ 



ds s'^P{s) - 1, (3) 



where (...) denotes frequency and ensemble averaging 
and A is the local mean eigenmode spacing. In case 
that propagation is not stopped by localization effects, 
var(s) « 0.273 for diffusive motion and var(s) = for bal- 
listic motion. By contrast var(s) = 1 indicates that a per- 
turbation cannot propagate indefinitely in the medium. 
A value of var(s) that increases (decreases) with system 
size signals that the perturbation will eventually (never) 
get localized fl^. 
The results for var(s) for the fBm and the Weiertrass 
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FIG. 3. var(s), as a function of the system size L, fixed disor- 
der a,, = 0.03 and u « 0.1 for a density profile described by 
a Gaussian white noise and a multifractal one generated by 
a multiplicative random cascade process (see text). In both 
cases the box dimension is D = 2 and the two point corre- 
lation function {Tj„rjo) oc 5{n). Since D > 3/2 all modes are 
localized. However localization effects are stronger (var(s) is 
larger) for the medium whose density is multifractal. 

density profile, at fixed w and disorder are shown in Fig. 
[Has a function of D. For D > 3/2 the variance var(s) in- 
creases with the system size. This indicates stop of wave 
propagation. By contrast for D < 3/2, var(s) decreases 
with the system size thus suggesting that, in this case, 
localization never occurs. Therefore in fractal media the 
transport properties for sufficiently weak disorder are to 
a great extent controlled by the fractal box dimension D 
of the density/stiffness. A natural question to ask is the 
degree of universality of this characterization in terms of 
D. In order to answer this question we investigate den- 
sity profiles with multifractal properties. In this case the 
box dimension can still be defined but the scaling proper- 
ties of the density profile depend on an infinite number of 
multifractal dimensions 2(f\. It is thus unclear whether 



other fractal dimensions control localization related ef- 
fects. 

Multifractal density generated by a cascade process.- We 
first study the effect of multifractality in a density profile 
whose two point correlation is {rjnrjo) oc 5{n). Our mo- 
tivation is to estimate the quantitative effect of a mul- 
tifractal media by comparing it with one, an uncorre- 
lated Gaussian noise, which is not multifractal but have 
the same two point correlation function. The multifrac- 
tal density fluctuations are generated iteratively by a 
variant of a multiplicative random cascade process [2l|. 
The length of the process doubles in each iteration. Ini- 
tially the density consists just of one value 77^°^ = 779. 
In the next generation it has two values, r;*^^) — m\ ijq 
and T^"' — 7TJ2 rjQ where nii are Gaussian random num- 
bers of zero mean and unit variance. After k itera- 
tions the density will have 2*^ different values given by. 
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?]i' '?^2;li ^'^'^ ^21' ~ '^i ' ^^21 ■ Following 



[22| we investigate the scaling properties of these data 
by using the multifractal detrended fluctuation analysis. 
The resulting density is characterized by a generalized 
Hurst exponent h(q) (see Q for a definition and Fig. 
I). For a monofractal curve, h((i) = H is independent 
of q. For this multifractal distribution, the box dimen- 
sion is D = 2 as that of a Gaussian uncorrelated disor- 
der. According to our previous characterization, local- 
ization should occur for any frequency and amount of 
disorder. A finite size scaling analysis (not shown), as 
the one performed in the monofractal case, confirms this 
prediction. This result might suggest that localization 
effects in both media are similar. However this is not 
the case. Multifractal corrections, which arises in higher 
order correlation functions, play an important role. In 
order to illustrate this we compare in Fig. [3] var(s) for 
both models as a function of L at fixed w and <Trj. We 
observe that: a) in both models var(s) increases with 
the system size. This is a signature of localized modes; 
b) Multifractal scaling enhances localization effects since 
var(s) is two order of magnitude larger in the this case. 
It is unclear to what extent this behaviour is universal. 
Multifractality leads to a broad range of scaling patterns. 
Therefore we cannot discard that in certain cases a sup- 
pression of localization effects can be observed. From our 
point of view the main conclusion of this analysis is that, 
in order to attain a good understanding of the transport 
properties in a highly inhomogeneous medium, is very 
important to know in detail the full distribution function 
that describes the medium density and not only its first 
moments. 

Multifractal density generated by a generalized random 
walk.- We now move to the case of a multifractal density 
in which h{2) can be tuned such that the box dimension 
D changes. Our motivation here is to test whether, as in 
the monofractal case, absence of localization is only ob- 
served for D < 3/2. For this purpose we study the mul- 
tifractal random walk introduced in [23]. The strength 
of disorder is set as in the monofractal case. Here we 
just state the definition of the model and refer to f23i| for 



mul 
'In 



i=l 



(4) 



where ry"^ is a fractional Gaussian noise, r{i) is a Gaussian 
correlated noise with (r^) — A^ In(Lo) and {r{i)r{j)) = 
A^ln/(|z — j\) with A a free parameter describing the 
strength of the multifractal scaling, Lq is the largest scale 
for which multifractal scaling is observed and f{\i — j\) = 
^ f.'^ . for \i — j\ < Lq — I and the unity otherwise. With 
these definitions the resulting fluctuations of the density 
profile jy"™' are multifractal with a set of multifractal 
dimensions that depends on A. As in the monofractal 
case we have carried out a finite size scaling analysis of 
the spectrum. The results for var(s) (see Fig. H]), indicate 
that propagation is arrested only for D > 3/2 since in this 
region var(s) decreases with the system size. Therefore 
our characterization for the existence of extended modes 
is still valid for media with a multifractal density. This 
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FIG. 4. var(s), as a function of D < 3/2 for uj ^ 0.1, 
(7^ — 0.1 and different L for a density profile described by 



a multifractal random walk with A 



0.1, Lo = 2'" and a 



monofractal with the same box dimension D. Multifractality 
enhaces localization effects but the critical box dimension is 
still D = Dc = 3/2. 

together with the previous results on monofractal media 
strongly suggests that the value of the box dimension D 
provides a characterization for the absence of Anderson 
localization in Id media. 

Finally we investigate whether our findings are consis- 
tent with the one parameter scaling theory [16]. Ac- 
cording to this theory var(s) can be expressed by a one- 
parameter scaling function var(s) = f{L/^(D)) where 
the scaling parameter ^ is the localization length in the 
localized regime, and the correlation length in the ex- 
tended regime. This relation implies that in a log-log 
plot of var(s) versus L all data should collapse in a com- 
mon curve when translated by an amount In^(D) along 
the horizontal axis. This curve has a single branch when 
there is no transition, while it develops two separate 
branches when a transition is present. In Fig. ([5]) we plot 



var(s) as a function of L for different values of D at a fixed 
w. The obtained result is fully consistent with the scal- 
ing hypothesis for var(s). Therefore the one-parameter 
scaling theory is still valid for this type of systems. 
Applications.- One of the main motivations to study 
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FIG. 5. (color online) var(s) a function of the system size for 
a; « 0.1 and cr^ = 0.1 for a density described by a fBm. Each 
colour/symbol corresponds to a different D. The observed 
data collapse is fully consistent with the one-parameter scal- 
ing assumption for var(s) [Ig| (see text). 



wave transport in fractal and multifractal media is its 
relevance in many realistic situations. We now address 
some of these applications. The distribution of fractures 
and densities in the earth inner structure is believed to 



be fractal/multifractal [ifl, [ill • Therefore a more com- 
prehensive understanding of the transport properties in 
such a medium would not only help characterize the earth 
inner structure but also would be relevant for practical 
problems such as the minimization of earthquakes effects 
and the detection of possible large-scale reservoir of oil 
and gas p4|. We note that our results provide a charac- 
terization of the conditions that a fractal medium must 
meet in order that a wave can propagate through it. 
A detailed understanding of the propagation of electro- 
magnetic waves in a highly inhomogeneous medium is 
a kev ingredient in the fabrication of different devices 
[3, Q (see introduction). In this case our characteriza- 
tion would be useful to optimize its capabilities. Likewise 
a comprehensive theory of the attenuation of acoustic 
waves in a fractal porous medium would improve dra- 
matically the design of sound absorbing materials. 
The study of matter waves in cold atoms settings could 
also benefit from these findings. In this case the system 
is a tight-binding Anderson model with a fractal-like po- 
tential. This type of potentials might be realized exper- 
imentally by using a ho log raphic mask combined with 
spatial light modulators 25 1. 

To conclude, we have studied wave propagation in a Id 
fractal media. Our main results are: a) the fractal di- 
mension of the potential D controls the strength of lo- 
calization effects; b) only for D < 3/2 wave propagation 
is not arrested by localization effects; c) multifractal cor- 
rections do not modify this result but have a profound 
impact in the transport properties, d) scaling theory still 
applies in these systems. 
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